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We introduce a computationally efficient approach to calculating the characteristics of excitons 
in quantum wells. In this approach we derive a system of self-consistent equations describing the 
motion of an electron-hole pair. The motion in the growth direction of the quantum well in this 
approach is separated from the in-plane motion, but each of them occurs in modified potentials 
found self-consistently. The approach is applied to shallow quantum wells, for which we obtained 
an analytical expression for the exciton binding energy and the ground state eigenfunction. Our 
results are in excellent agreement with standard variational calculations, while require reduced 
computational effort. 

PACS numbers: 71.35.Cc, 73.21. Fg, 78.67.De. 

I. INTRODUCTION 

Excitons play an important role in the band edge optical properties of low-dimensional semiconductor structures 
such as quantum wells, quantum wires, and quantum dotsi Quantum confinement of electrons and holes in such 
structures results in increased binding energy of excitons, their oscillator strength, and a life-time. As a result, 
excitons, for instance, in quantum wells, are observed even at room temperatures, and play, therefore, a crucial role 
in various optoelectronic applications To be able to calculate effectively and accurately exciton binding energies 
in the quantum heterostructures is an important problem, and therefore, a great deal of attention has been paid to 
it during several deca,desi 5 i 6 i 7 i 8 i 9 i 10 T 11 i 12 T 13 i 14 i 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 T 23 However, this problem is rather complicated, and 
while significant progress has been achieved, new material systems and new type of applications require more flexible, 
accurate and effective methods. 

Presently, the best results are usually obtained within the framework of the variational approach, where a certain 
form of the exciton wave function, depending on one or more variational parameters is being postulated. The exciton 
energy is then calculated by minimizing the respective energy functional with respect to the variational parameters. 
Unfortunately, even in the simplest (and, therefore, less accurate) realization of this approach it is not possible to 
express a value of the binding energy as a function of quantum well parameters. The best one can hope for is to 
obtain a set of complicated equations, which relate material parameters to several variational parameters. The latter 
are found numerically, and then used for numerical computation of the exciton energy. These difficulties result from 
the fact that variables in the Hamiltonian describing relative motion of the electron-hole pair cannot be separated: 
the presence of the quantum well potential breaks the translational invariance of the system, making it impossible to 
separate in-plane motion of the electrons and holes from the motion in the direction of confinement. 

The standard variational approach is based upon a choice of a specific form of a trial function, which is usually 
chosen in the form of a product of three termsi^LS The first two are one-particle one-dimensional electron and hole 
wave functions for confined motion across the quantum well. The third term describes the relative motion of an 
electron and a hole due to Coulomb interaction. The accuracy of the results depends on the complexity of the third 
term of the trial function and the number of variational parameters. The more parameters, the lower the binding 
energy, but, of course, the more extensive the calculations. There is ample literature dealing with accurate variational 
numerical calculations of exciton binding energy in quantum we llsi£i§iL2i2iiMLi2ii2iI4ii£ii&iLi2iiii. Most advanced of 
the calculations include the effects of Coulomb screening due to dielectric constant mismatch, as well as effective mass 
mismatch at hetero junctions and band degeneracyiSiiS 

Another approach discussed in the literature^2*2ii22iS is based upon an expansion of the electron-hole envelope 
wave function in terms of the complete system of eigenfunctions of a one-particle Hamiltonian describing motion of 
electrons (holes) in the respective quantum well confining potentials. The coefficients of this expansion represent the 
wave functions of the in-plane motion. They satisfy an infinite set of differential equations, which are coupled because 
of mixing of different electron and hole sub-bands induced by Coulomb interaction. Such a system can be solved only 
numerically after an appropriate truncation of the basis. Another way is to solve the system in diagonal approximation 
and to treat the off-diagonal elements with the help of perturbation theory. The improvement of accuracy in this 
approach faces difficulties related to the unknown errors due to basis truncation. 

Despite the significant achievement of the current approaches, they suffer from some principal limitations imposed by 
their very nature, and which cannot be, therefore, easily overcome. For instance, traditional variational approaches 
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are limited by the need to deal with a variational function of a particular form, which tremendously restricts the 
functional space over which the minimum of the energy is being searched. This problem cannot be circumvented 
by an increase in the number of the variational parameters because of the difficulties solving optimization problems 
with three or more parameters. Besides, calculations presented in most papers are not self-consistent (some limited 
attempts to introduce self-consistency, which were made in the pastiiiS are discussed below in Section II). At the 
same time calculations would become more important for material systems with wide band-gap materials. Thus, it 
is necessary to develop a method of calculating exciton binding energies, which would be more flexible and accurate 
than the existing methods, and which would allow to treat effects due to electron-hole interaction in a self-consistent 
way. 

In this work we suggest such an approach, which is based upon application of the ideas of the self-consistent Hartree 
method to the excitons in quantum wells. The idea of this approach is instead of imposing a particular functional 
dependence on the envelope wave function, to make a more flexible conjecture regarding the form of this function. 
Particularly, we present the total function as a combination of some unknown functions, which depend on fewer than 
the total number variables. Applying the variational principle to this combination we derive a system of equations 
describing both the motion of electrons and holes in the direction of confinement, and the relative two-dimensional 
in-plane motion of the exciton. Effective potentials entering these equations have to be found self-consistently along 
with the wave functions. 

This approach has a number of advantages compared to the previous methods. First of all, in its most general 
statement it must give better results for the exciton energy because we span a much larger functional space in the 
search for the minimum. Second, as it will be discussed below, this approach automatically gives a self-consistent 
description. Third, the approach naturally allows to incorporate external electric and magnetic fields, stress, disordered 
potential acting on electrons and holes in QW because of inherent inhomogeneities of structure. All these effects, 
which modify the single particle part of the Hamiltonian, appear automatically in self-consistent equations for the 
variational functions. We show that for the case of strong confinement in QW the self-consistent approach with 
factorized form of the envelope wave function allows one to achieve a very good agreement with the results of the 
variational method for the most elaborate trial functions used in the literature before. The relative motion of the 
exciton in the in-plane directions is described in this treatment by the Coulomb potential averaged with the wave 
functions of electron's (hole's) motion in the direction of confinement. The latter, in turn, is characterized by an 
effective confining potential, which is obtained by combining the initial quantum well potential with the appropriately 
averaged Coulomb potential. Unlike the perturbative method^, our approach takes into account the Coulomb mixing 
of the electron and hole sub-bands in a non-perturbative way, and is expected to give more accurate results even for 
the cases when such mixing is important. 

We show that the effective exciton potential has two different regimes of behavior. For large distances the potential 
has three-dimensional Coulomb tails, while at very small distances it becomes logarithmic as it would be for the true 
Coulomb potential of a point charge in two dimensions. A crossover between these two regimes takes place around 
distance r ~ d, which is the average electron-hole separation in the quantum well in the z-direction. 

Since the main purpose of this paper is to present the new approach and compare it with the results of calculations 
carried out by other methods, we chose a simplified model of a QW neglecting some of the effects that can be 
incorporated in the future. Moreover, we apply our approach to a particular case of a shallow quantum well, which 
allows for obtaining some analytical results, which give important qualitative insight into the properties of more 
generic models as well. We use a (^-functional model to describe shallow quantum wells that allows us to obtain an 
explicit form of the effective Coulomb potential. As a result, we derive a simple analytical formula for the exciton 
binding energy that depends only on one variational parameter. This formula gives results comparable to the best 
numerical results obtained by the standard variational approach. Thus we demonstrate that the method is both 
efficient and accurate; it can be applied to any quantum well with an average size of confinement smaller than the 
three-dimensional effective Bohr radius. 

The paper is organized as follows. In Sec. [HI we discuss the model and derive the general self-consistent equations. 
In Sec. Iflll we derive an analytical expression and discuss the different limits of the effective exciton potential. Section 
II VI presents the comparison of the binding energies results for the self-consistent approach and the standard variational 
method. The last section presents the conclusions of our work. The auxiliary details for the calculations can be found 
in three appendices. 

II. THE MODEL: 2D EXCITON IN SELF-CONSISTENT FIELD 

We assume that both conduction and valence bands are non-degenerate, and that they both have an isotropic 
parabolic dispersion characterized by the masses m e and mh (the heavy hole mass), respectively. Throughout the 
paper we use effective atomic units (a.u.), which means that all distances are measured in units of the effective Bohr 
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radius ob - 
where 1/fJ,* 



h 2 e//j,*e 2 , energies in units of n*e 4 /h 2 e 2 = 2 Ry, and masses in units of reduced electron-hole mass //*, 
= 1/m* + 1/to^. In this notation m e> /t = to* ^/a**, where to* ft are effective masses of an electron and a 
heavy hole. We assume that both the barrier and the well have close dielectric constants e as well as dispersion laws. 
Thus, we neglect a dielectric constant difference and an effective mass mismatch. One of the goals of this paper is to 
compare our method with existing approaches. Therefore, we have made some of these simplifications deliberately. 
Important effects such as valence-band mixing, non-parabolicity of the conduction band, dielectric constant and 
effective mass mismatches can be added to the model at later stages once the method is fully developed. 

After the standard procedure of excluding the center-of-mass of the perpendicular motion in the plane of the 
layersj£*£ the excitonic Hamiltonian is given by 
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where E g is a gap energy, 



z is the growth direction, r measures a relative electron-hole distance in the transverse 
direction r = y/ (x e — Xh) 2 + (y e — Uh) 2 , and Vi,2 are the quantum well confining potential in z direction for the 
electron and the hole, respectively. We have already assumed that the ground state must be independent of an angle 
in the xy plane, and excluded the corresponding term from the kinetic energy of the relative motion, K r . 

A variational principle can be used in two different ways for calculation of approximate solutions for the Schrodinger 
equation with the Hamiltonian |{T{|, The first approach is the standard variational method. It is well described in 
the literaturei£*§iL2ii2iiiii2ii2ii4ii£iI& According to this method, one needs to start from a variational principle for the 
functional E[W]: 



E[^f]= J y*HfdV 



with the additional normalization condition 



\WdV = 1. 



(2) 



(3) 



Then look for an approximate wave function within a class of functions of predetermined analytical coordinate 
dependence. These functions depend on several variational parameters, Ai, A2, • • •• Then the total energy 



E — E(Xi, . . . , A n ), 

and numerical values of variational parameters can be obtained from minimization conditions 

9£(A!,...,A„) 



0. 



1,2,. 



(4) 



(5) 



The success of the method depends essentially on the choice of the trial function. It must be simple enough to lend 
itself easily to the calculations, but must vary in a sufficiently large domain for the energy obtained to be closed to 
the exact one. 

Another way to calculate the approximate solutions of Eq. Q is to utilize the self- consistent approach. This 
approach also starts from the variational principle, Eqs. and However, instead of choosing a particular 

coordinate dependence of the trial function, we only assume a particular functional dependence on different coordinates 
for the entire wave function. Namely, we construct an approximate entire wave function ^(z e , Zh, r) with the help of 
the unknown functions ipi,ip2, ■■ ., where each function ipk depends on a lesser number of variables than the entire 
wave function. Considering variations of these functions independently, from the variational principle, Eqs. and 
P|). we obtain coupled integro-differential equations for ipk- 
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If localization in the quantum well is strong (the exciton "z-size" is smaller than its Bohr radius), then it is 
reasonable to suggest that the exact wave function for the ground state of Hamiltonian (|TJ is close to the simple 
product of functions of different coordinates 

^exactC?*, Z e , Z h ) ► ^trialO", Z e , Z h ) = lJj(r)Xe(Ze)Xh(z h ) . (6) 

Assuming normalization of every function in this product, we substitute function W in Eq. (J2J by the trial function 
ijBJl, vary each function in a product separately, and obtain the system of coupled integro-differential equations 

[K r + V r (r)] ^(r) = E X T/j(r), (7) 

[H e + V e (z e )] Xe(Ze) = E eX e{z e ), (8) 

[H h + V h {z h )]xh(z h ) = E h Xh(zh), (9) 
where the following notations for effective potentials are introduced: 

V r (r) = {XeXh\Vreh\XeXh), (10) 
V e ,h(Ze,h) = {lpXh,e\V re h\lpXh,e)- (H) 

The angle brackets imply that the integration of the Coulomb potential with corresponding wave functions is carried 
out over two of three independent variables. 

Solving system of equations JjJ— (JSJ we obtain the best approximation for the entire wave function in the form of a 
product © . The corresponding value of the total energy is given by Eq. that can be rewritten in the form 

E = = E e + E h + E X - (Xe\V e \Xe) - (Xh\V h \ Xh ). (12) 

The latter expression can be obtained by averaging each of Eqs. (J2J-© and adding them together. The electrostatic 
term between the electron and the hole is counted three times in the summations, and so has to be subtracted twice 
to give Eq. J5J). Thus the total energy is not just the sum of the exciton binding energy and the electron and the 
hole confining energies. The last two terms in Eq. (|12|l describe the renormalization of the total energy due to 
non-separability of the Hamiltonian. 

In order to solve Eqs. we a Pply the method of successive approximations. For strong localization inside the 

quantum well, the corrections to the single-particle energies due to effective Coulomb interaction potentials V et h(z et h) 

in Eqs. |(HJ| and JJjJ are small. Therefore, we begin by neglecting their contributions (^^=0) and solve the equations 

The obtained eigenfunctions are then substituted into Eq. (|1UII in order to get V r ° (r), a zero approximation for 
V r (r). This potential in turn should be substituted into Eq. (7J). The resulting equation, 

K r + V ( °\r)] V> (0) (r) = 4°V (0) W, (14) 

describes the formation of a two-dimensional exciton by an effective electron-hole interaction. The physical meaning of 
this effective interaction is a quantum mechanical average of the Coulomb potential with confinement wave functions. 
The ground state eigenfunction computed from Eq. (7J can then be substituted into Eq. to calculate a new 

approximation V ^ h {z e ,h) for the effective potentials. This process can be continued until the potentials are self- 
consistent to a high order of accuracy, i.e. until the condition 

MVVM » (Xe\V e \Xe) « (Xh\V h \ X h) (15) 

is fulfilled. Eqs. l(Tjl— l|12|l with the condition (|15() represent the complete system of equations for finding the minimum 
of the total energy for the Hamiltonian if the trial function has a particular functional dependence i|fj|. 

The described procedure has several advantages in comparison with the standard variational method. First of all, 
at each step we solve one-dimensional differential equations^ Second, even if the resulting Eq. (JJJ for the exciton in 
the effective field cannot be solved analytically, the explicit form of the effective potential i|10|) gives some additional 
understanding of the form of the exciton eigenfunction, and hence improves the accuracy of calculations. Finally, the 
convergence of the successive iterations itself allows us to estimate to what degree a given functional dependence of 
the trial function is close to the exact wave function. The energy difference between the successive approximations 
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shows how far the approximate energy is from the exact ground state energy. Certainly, we should expect a slow 
convergence for very broad and ultra-narrow quantum wells, where the entire wave function must be close to the wave 
function of the three-dimensional exciton. In this case, however, we can modify our self-consistent theory, rewriting 
Hamiltonian Q in terms of the new independent variables: the three-dimensional radius R (which is determined 
by R — \lr 2 + z 2 ), the angle 8 (that links z and R coordinates: z — Rcosd) and the coordinate Z of center-of- 
mass in the z direction [Z — (m e z e + rrihZh) /M]. Then we can apply the variational principle for a trial function 
Atrial #j Z) — ip(R)f(6)g(Z), and obtain corresponding self-consistent equations for functions i[>, f ,g. We do 
not want to dwell on this issue in this paper, since our primary interest is mainly focused on experiments where the 
average size of the particle localization in the z direction is smaller than the effective three-dimensional Bohr radius 
of the exciton. In the next section, we present the comparison of our self-consistent approach with the results of the 
standard variational method for the shallow quantum well. It is worth noting that the first iteration of our approach 
might be considered as an improved version of the standard variational method with a separable trial function © , in 
which functions Xe,h ar e the electron and the hole eigenfunctions without the interaction, and ip(r) is chosen in the 
form of a Is function for a two-dimensional exciton. 

We would like to note that some attempts to treat the Coulomb te rm in t he Hamiltonian JTJ in a self-consistent 
manner have been made in the pasti&iLi£*2ii22i2i For example, in Refs. ll7fl8l an incomplete self-consistent procedure 
for single-particle wave functions was perfor med. Due to the more complicated form of the trial function Atrial = 
tp(r, z e — Zh)Xe(z e )Xh(zh) the authors in Refs. [l7ll8l treated the first term i/j(r, z e — Zh) in the product by the standard 
variational method, adjusting variational parameters. Then an attempt was made to look for self-consistent corrections 
to the electron and hole wave functions, Xe,h, with the help of equations similar to our Eqs. iJSJ and J3J. In principle, 
it is possible to write the complete system of self-consistent equations for the trial function in the form ip(r, z e — 
Zh)Xe{z e )Xh{zh)- To do this, one needs to start again from the variational principle, Eqs. © and ©. Then by 
varying each function separately in the product one can get the complete system of the integro-differential equations 
similar to Eqs. 0-©- This system and the final expression for the total energy will have a more complicated form 
due to non-orthogonality of functions in the product for the entire trial function. 

On the other hand, the first iteration of our method results in Eq. I|14|) . which coincides with the zero approximation 
of the "truncated basis" approach^SLS when only the ground state confinement eigenfunctions are left. Our derivation 
of Eq. Q) shows, however, that it is more significant than merely a truncation of all but one term in the basis. The 
successive iterations of the Eqs. take into account the Coulomb mixing of the electron and hole sub-bands in 

a non-perturbative way, and give more accurate results even for the cases when such mixing is important. Numerical 
calculations confirm that, indeed, Eq. JTJ produces results, which are very close to those obtained by standard 
variational methods. 



III. EFFECTIVE POTENTIAL FOR EXCITON IN ^-FUNCTIONAL SHALLOW WELL 

We define a shallow quantum well as such a well, in which only one bound state exists for both electrons and holes. 
In general, the energy spectrum of the quantum well with height U and finite length L is E n — tt 2 x 2 1 / (2mL 2 ), where 
x n are the roots of the following transcendental equation: 

x = n — 2/tt arcsin ( x ) , n = 1, 2. . . . , iV (16) 

\Lu ) 

n — 1 < x n < n. 



Here we introduced a corresponding wave vector uq = \/2mU that characterizes potential height. The number of 
levels in the well is given by the condition 



N = l + 



u L 



(17) 



where [• • • ] denotes the integer part of the number. The condition given by Eq. I|17fl can also be interpreted as a 
condition of a new level appearing when the potential grows in the well. For example, the second level appears when 
Mo is equal to the wave vector of the ground state in the infinite well with the same length L: uq = n/L. 

A transcendental form of Eq. (|16fl as well as a piecewise character of the eigenfunctions present additional obstacles 
for further calculations of the exciton binding energies. Therefore, different approximations of the finite quantum well 
are often used. For a wide quantum well with several energy levels inside, the model of an infinite quantum well with 
a slightly larger effective length L e s = Ljx\ is an appropriate one^ It gives the same ground energy and correct wave 
function behavior. However, for shallow quantum wells with one level inside, the use of this model is not justified. 
Indeed, for an infinite quantum well the ground state energy grows with the decrease of the well's width, while the 
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FIG. 1: The energy levels dependence on the well's width. The ground and the first excited levels (solid lines) are shown for 
the finite quantum well. Dot-dashed lines show the ground state level for the infinite well of the same width and 5-functional 
potential I18H with a = UL. A vertical dashed line is the well's width at which the second level in the finite quantum well 
appears. It shows the range of applicability of the 5- functional potential with the effective strength a = UL c g. Parameters are 
taken for a conduction band electron in an Alo.3Gao.7As/GaAs quantum well (see below in a text). 



ground state energy in the finite well has a different dependence, and tends to the finite limit Uq when the width tends 
to zero: E\ — ► U — U 2 L 2 m/2. Narrow quantum wells have a different analytical limit of the ^-functional potential^ 



V(z) = U-a5(z), (18) 
where a is a ^-potential strength. If we define this parameter as 

a = UL cS =J— (19) 

where L G ff is chosen to match the ground state energy of the finite well problem, then the well-width range of 
applicability of this approximation is extended up to the moment of the appearance of the second level in the finite 
quantum well. For typical parameters in AlGaAs/GaAs structures it corresponds to a well size L « 40 A. Obviously, 
L e fi — > L when L tends to zero. Figure ^ shows typical energy dependence on well's width for the electron in the 
AlGaAs/GaAs quantum well for the finite quantum well and its approximations. Comparing the behavior of curves 
for the ground state of the finite width well and its 5-functional approximation with strength a — UL, one can see 
that ^-functional curve stays always on the left. It means that the effective length parameter, determined by Eq. i|19|) 
should be smaller than the actual well width, which is opposite to the case of effective infinite quantum well width. In 
some sense, the model of 6- functional QW is complimentary to the model of effective infinite quantum well24 (EIQW), 
which is used to approximate finite QW with large widths (and/or barrier heights), when the number of levels in a 
well is large. Indeed, the more discrete levels exist in the QW the better the EIQW model works, but it fails gives 
a wrong eigenstate dependence on L, when the well has only one level. On the other hand the 5-functional QW is 
not applicable for quantum wells with more than one level. The delta-functional approximation is applicable both to 
very narrow QW and to wells with a small band-gap offset^ i.e. when the well width and/or the band offsets are 
very small so that the carrier wave functions are mostly in the barrier region. For the 5-functional potential it is more 
convenient to count energies from the barrier band edge rather than from the bottom of the well. In terms of the total 
Hamiltonian it means that the energy band gap constant is the barrier's energy band gap: Eg = Eg 1 + U e + Uh- 
The energy and wave function of a single localized state are well-known: 

Ee.h = ~7T^ ' Xe.h(z) = y/KTh 6Xp(- K e h\z\) . (20) 

Parameters K e h = Oi e h m e,h determine the localization of wave functions of an electron and a hole, respectively. It 
is worth noting that even for a very shallow quantum well, the localization length (~ 1/k) might be much less than 
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FIG. 2: An effective self-consistent potential profile for a two-dimensional exciton. The solid thick line represents Eq. 
(121 1 . The dashed-dotted line is the approximation (r 2 + d 2 ) _1//2 , the dashed line is the logarithmic regime Eq. 12511 for small 
r, the dotted line represents 1/r behavior. All data are for d/as = 0.16 corresponding to L = 20A finite quantum well in 
AlGaAs/GaAs materials. 



the effective Bohr radius. In this case, we can expect a quasi-two dimensional behavior for the exciton, justifying the 
approximation for the mean field function in form In the case of the AlGaAs/GaAs quantum wells, the electron 
(hole) localization length is smaller than Bohr's radius up to L 5A. 
With the help of the wave functions (|20|l . it is possible to obtain the analytical expression for the effective exciton 

potential V r (Y). The details of these calculations are given in Appendix lAl The result is 

V { r °\r; K e ,Kh) = — 2 2 KeKh 2 [n h T(n e r) - n e T(n h r)] , (21) 

K h K e 



where the function T(nr) is a combination of zeroth-order Struve and Neuman functions 
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T{nr) = ~ [H (2/tr) - Y (2/tr)] . (22) 

The behavior of the potential (|21|l is shown in Figure [3 

The behavior of the potential l|21|) has two regimes that are determined by the parameter 









li 


We ' 





(23) 

This parameter has the meaning of an average electron-hole separation in the z direction. For large distances, potential 
Q21JI has asymptotic behavior 

V4 sym (r) « -^=L=, r>d. (24) 
At small distances, the attraction becomes stronger. It has logarithmic behavior: 28 

V sm (r) w \ [Hr/d) - 1 + 7] , r < d, (25) 
d 

where 7 = 0.5772. 

Thus, the effective electron- hole interaction for the exciton in the quantum well starts from the true logarithmic 
Coulomb potential of a point charge in two dimensions that smoothly transforms at distances r ~ d to the screening 
potential l|24|l with three-dimensional Coulomb tails. For the strong confinement d< 1 we can approximate V e s(r) 
by Eq. (|24H for all distances and take into account the logarithmic part on the next step as a perturbation. 
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It is interesting to note that potential 1(24(1 can be obtained without the self-consistent procedure from the following 
simple intuitive consideration. At the first step, lets us neglect the electron- hole interaction in Hamiltonian (|T|> . Then, 
we can solve the one-dimensional one-particle Schrodinger equations in the quantum well, and find the average square 
of the distance between the electron and the hole as 

d 2 = (Xe{z e )Xh(zh)\(z e ~ Z h f \ Xe (z e )Xh (z h )) ■ (26) 

This yields the same result as Eq. (123(1 . The next step in the approximation of the Hamiltonian is to include the 
Coulomb attraction term, where (z e — z^) 2 is substituted by its average value ( (z e — z^) 2 ). 



IV. NUMERICAL RESULTS: COMPARISON WITH STANDARD VARIATIONAL APPROACH 

The Schrodinger equation for a radial wave function in the central field 1(21|) does not have an analytical solution. 
To obtain an approximate analytical expression for the exciton binding energy we will use the following procedure. 

First of all, let us stress again that despite the fact that we consider a shallow quantum well with one single particle 
eigenvalue inside, the strong confinement persists up to very small widths. The parameter d is a good indicator of 
such confinement. For example, in Table[l]the data are presented for Alo.3Gao.7As-GaAs materials. We can see that 
for the well's width of lOA this parameter is about one quarter of the three-dimensional Bohr radius and even smaller 
for larger quantum wells. For the case of strong confinement, d <C f , the effective potential l(21|l can be represented 
by Eq. 1(24(1 almost everywhere. Therefore, at the first step, it is reasonable to substitute the potential 1(21(1 by its 
asymptotic form 1(24(1 for all distances. 

Appendix 151 yields the details of numerical calculations and analytical limits for the ground state in the potential. 
Although the Schrodinger equation with the potential 1(24(1 also does not have an analytical solution, we discovered 
that the ground state energy, obtained by the variational method for the single parameter trial function 



Atrial 



2exp(d/A) 

y/\(\ + 2d) 



exp(- 



<P/\), 



(27) 



coincides with the exact one with excellent accuracy. To check this, we performed a precise numerical integration 
of the Schrodinger equation based on Pruefer transformation and a shooting method. The difference in the ground 
state energies for the whole range of the parameter d was less than 0.01%, or ~ 10~ 3 meV! Such an agreement can be 
explained by the fact that trial function ((27(1 has a correct analytical behavior for both small and large distances, r. 
The expression for the ground state energy obtained for the trial function 1(27(1 is given by 



A 1 + 2d/\ 



1 

A 2 



(2rf/A) 2 £i(l, 2d/\) exp(2d/A) 
1 + 2d/\ 



(28) 



where Ei(x) is the exponential integraliS 5 . The variational parameter A changes from 1.1 to 0.74 when the average 
electron-hole distance d varies from 0.11 to 0.48. The latter corresponds to the quantum well widths range from 40A 
to 5A for the AlGaAs/GaAs structures. The behavior of the parameter A as a function of d is shown in the insert of 
Fig. |3 At small d it has the following form: 



A 



4d- 



A 2 



4 In 



12A + 4 7 + 1 



(29) 



At small distances the effective potential 1(21(1 differs from Eq. ((24(1 . The correction to the energy due to this difference 
can be taken into account with the help of perturbation theory: 



E 



(ip trial \AV\ Atrial) 



2d 
A 2 



1 2d {2d 
-- 0.557— + 0.563 — 

2 A \ A 



i 



\n{r/d) 



■ dr 



(30) 
(31) 



The last column of the Table |U represents the final sum E^ = E 



(o) 
J x 



+ 2#> 



for the exciton binding energy obtained 



by the self-consistent approach. To check the accuracy of our method we compared the results of the self-consistent 
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FIG. 3: The exciton binding energies in the shallow quantum well for different well widths L. The solid line represents the 
two-dimensional exciton ground state, Eqs. 1281 and 1301 in the effective potential. The dotted-dashed line is the exciton 
binding energy obtained with the help of the variational method with the trial function 132H . The insert shows the dependence 
of the parameter A on the average electron-hole separation d in Eq. 1281 . Both parameters are expressed in Bohr radius units 
and are presented for the same range of quantum well widths. 



TABLE I: Effective parameters of single quantum well potentials, electron, hole and exciton binding energies for different 
quantum well widths, L. Single particle electron and hole energies, E®\, are the initial step in the self-consistent field iterations 
[see Eq. 1131 1 when effective Coulomb terms V e ,h are omitted. These energies determine the strength parameters, L^, for 
corresponding 5-potentials and the average distance, d, between an electron and a hole [Eq. 1231 1. Energies ^W^ 2 )^ 3 ' are the 
exciton binding energies obtained by variational method with the trial functions given by Eqs. 13211 . 1331 . and 1341 . respectively. 
The last column shows the energy Ex = Ex + E^ [see Eqs. 1281 . 131)1 1 obtained by the first iteration of the self-consistent 
approach. The calculations are based on the following physical constants* 2 — m* = 0.067mo, m* h — 0.45mo, U e = 340 meV, 
Uh = 70 meV, e = 13.8. For these parameters the effective Bohr radius is as = 1254 and the energy atomic unit is equal to 
8.33 meV. 



L(A) 


icff (A) 


(A) 


d(a.u.) 


-Ei 0) (a.u.) 


-4 0) (-U.) 


-E^ (a.u.) 


-Ef (a.u.) 


-Ef (a.u.) 


-E% (a.u.) 


10 


9.15 


8.89 


0.257 


5.11 


1.37 


1.0799 


1.0975 


1.1041 


1.0948 


20 


15.16 


14.12 


0.158 


14.02 


3.46 


1.2617 


1.2779 


1.2843 


1.2754 


30 


18.63 


16.88 


0.130 


21.18 


4.95 


1.3314 


1.3467 


1.3527 


1.3443 


40 


20.70 


18.43 


0.118 


26.12 


5.90 


1.3655 


1.3802 


1.3855 


1.3779 



approach with the results of the standard variational method for three different trial functions. These trial functions 
have the following forms: 

ip^(z e ,z h ,r; A) = ^K e K h exp(-K e |z e | - K h \z h \) exp(-r/A), (32) 
ipW(z e ,Zh,r;\,0) = y/K e K h exp(-K e |z e | - nh\zh\) exp(- \J r 2 +/3 2 /A), (33) 
^ (3) (z e ,^,r; A,/3) = ^n e K h exp(-K e |z e | - n h \z h \) exp(- yjr 2 + /3 2 (z e - z h ) 2 /X), (34) 

The first two functions are separable, while the third one is non-separable. The first wave function has one variational 
parameter A, and two others have two variational parameters A and (3. 

The details of variational calculations are given in Appendix[n] Our results and their comparison with data obtained 
by the standard variational procedure are presented in Table [i] 

We can see that the single parameter separable trial function, (z e , Zh, r; A), gives higher binding energy with 
maximal relative difference of 1.3%. The non-separable trial functions have slightly lower binding energies than the 
first iteration of the self-consistent approach. The maximum relative difference between the self-consistent binding 

(2) 

energy Ej£ and the binding energy E x is 0.24% (for AlGaAs/GaAs structures it corresponds to a difference of 0.02 
meV). For it is 0.84% (0.07 meV). These results demonstrate that even the first iteration of our approach gives 
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an excellent agreement with the variational results. Subsequent iterations further decrease the ground state energy. 
Preliminary research on the convergence of the self-consistent approach shows that the next iteration gives the value 
of the binding energy lower than the variational calculations with trial functions 1)33(1 and 1)34(1 . These results will 
be published elsewhere. However, in the particular case considered here the difference by 0.07 meV between more 
complicated variational approach and the results presented here is already much smaller than the accuracy of, for 
instance, optical absorption experiments^, which is of the order of 2 meV. From this point of view the discrepancy 
between the two methods is negligible. At the same time, our approach is about 100 times faster than the variational 
method even with the simple enough trial function, Eq. I(33|) . It also helps to avoid a numerically difficult task of 
finding minima of several non-polynomial functions. The additional physical information about the effective potential 
allows one to find one of the variational parameters of function 1(33(1 without minimization. 



We introduce a self-consistent approach for calculations of the exciton binding energy in a quantum well. For the 
case of strong confinement, the self-consistent Hamiltonian is separable and consists of three parts: one-dimensional 
Hamiltonians for electron and hole confined motions across the quantum well, and the Hamiltonian, which describes 
the motion of a two-dimensional exciton in the effective central field potential. This effective potential is a result of 
averaging over z coordinates of the Coulomb interaction and the quantum well potential. As a function of distance 
the effective potential has two different regimes of behavior which are determined by the average distance between 
electron and hole inside the quantum well, d. For small distances, r < d, the effective potential has a logarithmic form 
of a Coulomb potential of a point charge in two dimensions. At a distance r ~ d this behavior crosses over to the 
three-dimensional Coulomb screened potential, — 1/ \fr 2 + d 2 . For the case of the shallow quantum well, analytical 
formulae, Eqs. 1(28(1 and 1(30(1 . for the exciton binding energy are obtained. Even though the use of these formulae 
requires computational time which is by orders of magnitude smaller than that of standard variational calculations, 
the results obtained by both methods are in an excellent agreement. The differences between the exciton binding 
energies in two approaches are generally smaller than 1%. For AlGaAs/GaAs structures it corresponds to differences 
smaller than 0.1 meV, i.e. smaller than corrections^ due to non-parabolicity of the bulk conduction band or dielectric 
constant and effective mass mismatches. One can expect that the developed method can lay the foundation for models 
that incorporate these important effects as well as take into account an influence of the external electric field. 

All results discussed above are obtained for the first successive iteration. Obviously, the next iterations will lower 
the ground state energy even more. For example, the next step is the substitution of the wave function 1(27(1 into 
Eq. ljTT|) . which is approximately equivalent to the appearance of additional oscillatory potentials in Eqs. (JSJ and J*jJ: 



These potentials will slightly change the single-particle energies and will localize the tails of the eigenfunctions of the 
electron and the hole ground states. The latter will be manifested in a decrease of the average electron- hole distance 
d, and, therefore, will lower the exciton binding energy in the next successive iteration of Eq. l|7jl. 

In derivation of our results we used three different approximations. They are: (i) the self-consistent approach itself, 
(«) the use of the factorized form of the wave function, (iii) and the use of the delta-functional potential for a shallow 
quantum well. The self-consistent approach is broader than the standard variational method since we do not have to 
specify a particular functional dependence for a trial function. Instead, we just suggest that the trial function consists 
of some combination of unknown functions. The factorized form is the simplest form for such a combination but it 
is not required by our method. The method can be applied to other physical models where different types of trial 
functions would be more natural. For example, for a wide double quantum well, the better choice of the trial function 
of the ground state is a superposition of two factorized single-well functions. It will result in a system of coupled 
equations similar (but more complicated) to Eqs. Q 0. This approach can also be straightforwardly expanded to 
include other effects, which were neglected in this paper. For instance, in order to take into account the dielectric 
mismatchj2ii2iii we would need to correct the expression for the effective potentials, Eq. ((10(1 and Eq. (|1 1|) . including 
effects of image charges into the respective integrals. The valence band degeneracy and anisotropy can be included 
by introducing a four component trial function, for which the self-consistent equations will have the similar form as 
Eqs. Q-01, but they should be understood as matrix equations. 

Turning to the particular factorization used in our problem, it is obvious that the more strongly the exciton is 
localized in the z-direction (size of wave function in z-direction compared to the three-dimensional Bohr radius) 
the better our approximation works. If we consider the case of Alo.3Gao.7As-GaAs quantum well, it means that 
our method will work for any quantum well with a width less that 100 A, but the best convergence will happen 
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(35) 
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somewhere around 30 — 50 A. Correspondingly, the factorized form of the trial function for the self-consistent method 
is applicable for any structure (e.g., asymmetric quantum well or quantum well in electric held), if a wave function 
localized in z-direction has such an extension. In a similar matter the self-consistent approach can be applied to the 
lower dimension systems such as quantum wires and quantum dots. Moreover, preliminary consideration showed that 
the problem of divergency of the exciton ground state in a one-dimensional Coulomb potential, which arises in other 
approaches to quantum wires, in this method does not appear at all. 

The J-functional potential is a good approximation if a one-dimensional quantum well has only one level (is shallow) . 
For a typical case of Alo.3Gao.7As-GaAs quantum well, it gives an applicability range of L < 40 A. The approximation 
of the 5-functional potential gives simple single-particle wave functions that significantly simplify calculations of the 
effective potentials in Eqs. IjlOfl and (|lll) . If asymmetric quantum well or double quantum wells, or quantum well in 
electric held are shallow, the delta-functional approach will be applicable to such models. An asymmetric quantum 
well can be modelled by a quantum barrier and the ^-functional potential (S(z) — ► <5(z)[l + A6(z)], where 9 is the step 
function). An advantage of using the (5-functional potential is especially clear in the case of the quantum confined 
Stark effect for a shallow quantum well, where it allows one to obtain additional results related to the held-induced 
resonance widths for a single-particle as well as for the exciton quasi-bound states. These results will be published 
elsewhere. 
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APPENDIX A: EFFECTIVE POTENTIAL FOR 2D EXCITON 



The effective held for quasi 2D exciton is given by the integral 



/OO r-C 
dz e I 
-OO J — I 



dz h 



We{Ze)\ 2 Wh{zh)Y 



yjr 2 + (z e - Z h ) 2 

For a shallow quantum well approximated by the ^-function potential it gives 

2kiK2 2ki«2 

■s/r 2 + (Z! - z 2 ) 2 yjr 2 + {zi + z 2 ) 2 
After making the coordinate transformation £ = z\ — Z2, rj = Z\ + z 2 and taking into account that 



V eS (r) 



dz\ \ dz 2 e 



= V X + V 2 . 



POO POO -i POG pTj 1 POO pOO 

/ d Zl dz 2 f(t, v ) = - dr, d4m,r,) = - dU dn [/(-£, v) + mv)]> 

Jo Jo 1 Jo J-n z Jo J£ 



(Al) 



(A2) 



(A3) 



the hrst integration in these two interated integrals becomes trivial and the second integration can be expressed 
through the function 



T(nr) 



cxjp(—2nt)dt 7T r , 
%= = = ~n [Ho(2 K r) - Y (2 K r)} 
o Vr + t z i 



(A4) 



where Ho is the zeroth-order Struve function and Yo is the zero-order Neumann or Bessel function of the second 
kindi2& Then potentials V\ and V 2 can be expressed as 

KiK 2 



and the hnal result yields Eq. I|21() 



V c s(r;ni,K 2 ) 



Kl ± K 2 

1k,\k 2 



[T( K2 r) ± T{ Kl r)} 



K 2 K l 

In the case when K\ = k 2 = k this expression is reduced to 

Ves(r; k) = -K [T(nr) - nT'{nr)] , 

where 

T'(Kr) = ^ = 2r [^Y 1 (2 K r) - jHx(2 K r) + 1 
ok L 2 2 



(A5) 

(A6) 

(A7) 
(A8) 
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d/a B r/a B 

FIG. 4: (a) The ground state energy for the potential — (r 2 + d 2 ) _1//2 , Eq. 1B3L as functions of the parameter d. Solid line, 
representing numerical results, and dashed line, obtained by the use of variational formula 1B5II . are practically indistinguishable. 
The dot-dashed lines correspond to perturbation theory for small d and an asymptotic oscillatory energy — eP 1 + d" 3 ^ 2 for 
d ^> 1. (b) The ground state wave function for d = 0.2. The numerical results are shown by the solid line. The dashed line 
corresponds to Eq. lB4t for normalized trial function. 



APPENDIX B: EIGENVALUES AND EIGENFUNCTIONS IN 2D CENTRAL POTENTIALS 



The Schrodinger equation for a radial wave function in a central field is 



d 2 1 d m 2 

+ T -2U(r) + 2E 

or z r or r z 



ip(r) = 0. 



For the Coulomb potential, U(r) = —Z/r, the ground state is well-known: 



E x = -2Z 2 , <p x {r) 



AZ 



-2rZ 



2tt 



Let us consider now the potentia 



26 



U(r) = -l/y/r*+<P. 



(Bl) 



(B2) 



(B3) 



The ground state energy and the eigenfunction for Eq. (|B1|) with the potential (|B3|I are presented in Fig. 

It turns out that the numerical results are in an extremely good agreement (~ 0.01%) with the results of variational 
method for the trial function 



2exp(d/X) , ^,, s 

Lptrial = 7fn WT\ exp ("V r + d / A )' 

y \\\ + 2a) 



(B4) 



where A is the only variational parameter. The corresponding value for ground state energy is obtained from mini- 
mization of the functional 



E X (X) = ~ ' 



1 



A 1 + 2d/\ A 2 



1 - 



(2d/A) 2 £i(l, 2d/X) exp(2d/A) 



l + 2d/A 



(B5) 



Such a good agreement can be explained by the fact that the trial function i|B4|) has correct asymptotic behavior in 
both limits of small and large r. 

Far away from the "atomic residue" U(r) ~ —1/r, the wave function for the s-state obeys the Schrodinger equation 

A 1 a 2 

- 1 v--v = - T v 

with the solution 



(r) = A a r^ a -^ 2 e- ar , 



(B6) 
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where terms of the order of tp/r 2 are neglected. Here A a and a are atomic parameters. Their magnitudes are 
determined by the behavior of an electron inside the "atom" . 

Simple analytical estimates for the cases of small and large d for ground state energy and the asymptotic coefficient 
A a can be made. When d -C 1 the potential is only slightly different from the Coulomb potential and the contribution 
to the energy can be obtained using perturbation theory: 

E x = -2+ (0|V|0) w -2 + 16d(l + 2d\n(2d)) . (B7) 



Note, that due to the big factor in front of d in Eq. I|B7[I. the convergence radius for the perturbation series is small 
enough. 

Assuming that the asymptotic form (|B6|) of the wave function is valid for all values of r, we obtain the following 
simple estimate for A a 

roo 2 1 /a-l/2 l/a+l 

2«Al / r 2 ' a exp(-2ar) dr = 1 =>- A a = . (B8) 

In the opposite case, when d^S> 1 the solution has to be close to the oscillatory one 

E x = - l - + ^ I - 2 {n + l), (B9) 

where n = 0, 1, 2.... The main contribution to the normalization comes from the Gaussian part of the wave function 
because the contribution of the Coulomb tails is negligible. One obtains for the state with n — 



where w = 1/d 3 ^ 2 . Then, 



tp(r) « <Posc(r) = J - exp (-ujr 2 /2) , (BIO) 



A * ^ (B11) 

APPENDIX C: VARIATIONAL METHOD FOR EXCITONS IN 5-FUNCTION QUANTUM WELLS 

Following standard procedures the envelope variational exciton wave function in a quantum well can be presented 
as a product of three terms, 

*(z e) Zh, r; A,) = Xe{z e )xh(z h )(j){r, z ei z h ; A,-), (CI) 

where Xi is the set of variational parameters and <f> is the variational wave function which minimizes the total energy 
of the Hamiltonian Q. Two other factors Xe,h{Ze,h) are simply normalized eigenfunctions of the one-particle electron 
or hole Hamiltonians of the quantum well: 



E e ,hXe,h(z e ,h) 



1 d 2 

2 ate t h5(Ze,h) 



Xe,h{Ze,h), 



2m e ,A dz 2 e h 

Xe,h( Z e,h) = y/Ke,h exp (- K e ,ft | Z e ,h | ) , 

K e ,h = m eih a eth , E eJl = K 2 eh /2m eJl . (C2) 
To obtain more confident results we made calculations with three different trial functions ip: 

V; (1) (r;A) = exp(-r/A), (C3) 



^ (2) (r;A,/3) = exp(- ^r 2 + f3 2 /X), (C4) 
^ 3 \r,\z e ~z h \;X,(3) = eM-^/r 2 + f3 2 (z e ~ z h ) 2 / X), (C5) 

The first two functions are independent of z coordinates, while the third one is non-separable with respect to z. The 
first wave function has one variational parameter A, and two others have two variational parameters A and [3. 



14 



The total exciton energy is the minimum of the functional 



E = 



(C6) 



For the first two trial functions, which are independent of z, the functional l|C6(l can be further simplified. In this 
case, the energy can be presented as 



E = E e + E h + K/N + V/N, 



where 



N = < 



*|*),= (V(r)|VW) 
1 
2 

1 



K = (${r)\K r \il>(r)) 
V = 



^Jr 2 + (z e - z h ) 2 

Calculations for the trial function ip^- 2 ' yield 

(a) = exp(-2/?/A)A(A + 2/3) 
4 



dA, 

|*) = f (V cS (rM(r)\ 2 ) dA, 



-(2) = exp(-2/3/A) 



1 + 



2(3 (2(3 



A V A 



E 1 (2/3/A)exp(2/3/A) 



v (2) = _^1^2_ [K2y(Ki5A;/3) _ Kiy(K2jAj/3)]5 



where Ei (x) is the exponential integral^ and 

r-rr/2 



V(k,X,/3)= cos(</>) d4> / exp(-2/ci?sin(0) - 2^-R 2 cos 2 (0) + /3 2 /A)i? eLR. 

J o Jo 

For the first trial function ip^ 1 ' (r; A) = ip^ (r; A, 0) all integrals have analytical expressions: 
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1 \/1 + A 2 k 2 + 1 
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(C7) 

(C8) 
(C9) 

(CIO) 



(Cll) 



(C12) 



(C13) 



Integrals for the non-separable trial function can be numerically estimated following the procedure described in 
Ref. lllH 
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